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1. Introduction 

The majority of the stellar-mass black holes (BHs) present in the modern universe are expected 
to form by the collapse of protoneutron stars in core-collapse supernovae (CCSNe). There are 
several channels through which this can occur: ( i) the supernova shock fails to be reenergized and 
accretion pushes the PNS over its maximum mass, leading to collapse (c.f. [11, 10, 3, 6]), (ii) 
fallback after a successful CCSN can lead to PNS collapse [13], and (Hi) nuclear phase transitions 
or PNS cooling may cause the PNS to collapse during its cooling phase (e.g., [1]). 

In this contribution to the proceedings of the 11th Symposium on Nuclei in the Cosmos, we 
investigate several aspects of BH formation in CCSNe that ai^e related to thermal effects. In §2, 
we briefly outline the simulation code we use to study core collapse and BH formation and dis- 
cuss our implementation of nucleon-nucleon Bremsstrahlung, a recent addition to our code. In 
§3.1, we show that in failing CCSNe this interaction increases the neutrino luminosity, removes 
thermal support, and causes the PNS to collapse to a BH earlier than when the interaction is not 
included. Finally, in §3.2, we artificially drive explosions in which the resultant PNS is thermally 
supported above the maximum baryonic PNS mass. Using this approach with an extreme model, 
we demonstrate delayed BH formation due to PNS cooUng. 

2. Methods 

2.1 GRID 

We use the open-source, spherically-symmetric, general-relativistic (GR) Eulerian hydrody- 
namics code GRID, which is available at http: / /stellarcollapse . org. We refer the 
reader to Ref. [5] for the full details of GRID and give here only a brief description. GRID follows 
the formalism of Romero et al. [7], which uses the radial gauge and polar slicing gauge condition 
in spherically-symmetric GR. This reduces the metric to a diagonal form with two metric functions 
that are determined by hydrodynamic variables. GRID is designed to follow the evolution of stars 
beginning from the onset of core collapse to BH formation. Its GR hydrodynamics module utilizes 
high resolution shock capturing methods. GRID makes use of several microphysical equations of 
state (EOS), also available online at http : //stellarcollapse . org, and uses a 3-neutrino 
leakage/heating scheme (described in §2.2). Here, we make use of the EOS from Lattimer & 
Swesty (I99I) [4] with an incompressibiUty ofKo = 220 MeV (LS220). 

2.2 Neutrino Treatment & Nucleon-Nucleon Bremsstrahlung 

GRID uses a hybrid neutrino leakage scheme based on the schemes of [9] and [8] and includes 
the standard reactions therein. It includes three neutrino types, Ve, Ve, and Vx = {v^, v^,Vt,Vt}. 
The neutrino emission rates, both energy and number, are determined from an interpolation be- 
tween the diffusion and free streaming limits. These emission rates are combined with a param- 
eterized neutrino heating scheme to calculate the total energy emitted from a zone. The parame- 
terization of our heating scheme, via the scale factor /heat> allows us to drive artificial explosions. 
Full details of GRlD's neutrino leakage/heating scheme can be found in [5]. An application of our 
parameterized heating scheme can be found in [6]. 
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A process for neutrino pair production not included in theleakage schemes of [9, 8] is nucleon- 
nucleon Bremsstrahlung (NNBrems), i.e. N + N N + N + vv. NNBrems dominates the neu- 
trino pair production at high matter densities. In CCSNe, this regime is very optically thick and 
NNBrems is a subdominant contributor to the total neutrino luminostiy. However, during PNS 
cooling or NS-NS or NS-BH mergers, emission from NNBrems may be considerable. As a ther- 
mal neutrino pair-production process, we include NNBrems into GRID to investigate its effect on 
BH formation. We use the single neutrino pair emissivity for NNBrems from [2] ^ , 

e„, = 2.0778 xl03°c(^,?+^p + y^«^p)pu(^) ergscm-^s-^, (2.1) 

where pi4 is the matter density in units of 10^"^ g cm^^, T is the matter temperature in MeV, X,i 
and Xp are the mass fractions of neutrons and protons respectively and ^ is a factor that absorbs all 
interaction ambiguities. We adopt the value of ^ = 0.5 following [2]. For GRID, we require the 
single-neutrino emissivities (both energy and number). Absorptive opacities are not needed (but 
could be obtained via Kirchhoff's law), since GRlD's leakage scheme includes opacity contribu- 
tions only from scattering and from chai^ged-current absorption processes. With the publication of 
this proceedings article, we have updated GRID (to vl.03) to include NNBrems as an option in the 
leakage scheme. 

3. Results 

3.1 Nucleon-Nucleon Bremsstrahlung 

We perform 12 simulations in which we investigate the effect of NNBrems in BH-forming core 
collapse simulations. We choose 6 progenitor models from the "w series" of Woosley et al. (2002) 
[12], namely the 35, 40, 45, 50, 60, and 75M0 models with metalhcities of IO^'^Zq. These models 
have large iron cores at the onset of core collapse and, as we shall see, form BHs very soon after 
bounce. They yield large temperatures in their PNSs and therefore are ideal for studying thermal 
effects. For each model we run simulations with and without the NNBrems emissivities (Eq. 2. 1) 
included in our neutrino leakage scheme. Table 1 presents the results of these simulations, we 
include the time to BH formation and the total neutrino energy emitted between bounce and BH 
formation. 



Model 


u35 


m40 


m45 


m50 


m60 


ulS 




^BH [s] 


0.522 


0.478 


0.578 


0.747 


0.630 


0.311 


NNBrems 














V- 


energy [B] 


257 


260 


111 


281 


233 


184 


No NNBrems 

V- 


?BH [s] 
energy [B] 


0.526 
249 


0.482 
253 


0.582 
268 


0.749 
273 


0.635 
229 


0.311 
179 



Table 1: Effect of NNBrems on time to BH formation and total emitted neutrino luminosity until BH 
formation. The progenitor models are taken from the lO^^Z© model set of [12] and simulated with the 
LS220EOS. 1B= lO^'ergs. 

'[2] mistakenly included the factor in the numerical constant as well as in the expression for Q^i,. We use the 
corrected expression here (Burrows, pnva?e communication, 2010). 



3 



Thermal Effects on Black Hole Formation in Failed Core-Collapse Supernovae 



Evan O'Connor 




Figure 1: Shock radii (left panel) for the m75 model of [ 12] with 7 values of /heat- For each value of /heat, we 
plot in the right panel the baryonic mass of the PNS after accretion has ceased versus the time each model 
remains stable. Less massive PNS remain stable for longer times. The dotted line indicates the maximum 
cold neutron star baryonic mass for the LS220 EOS. 



NNBrems results in a slightly higher^ (~3%) average neutrino luminosity as is demonstrated 
by the higher total energy emitted in neutrinos. This is due to the increased neutrino emissivity in 
the dense core of the PNS where NNBrems dominates. The extra loss of thermal energy from this 
region reduces the thermal pressure, thus removes thermal pressure support from the PNS. This, 
in turn, lowers the maximum PNS mass and causes the dynamical collapse of the PNS to begin at 
somewhat earlier times (At ~ 0(ms)) compared to the case with no NNBrems. 



3.2 Thermal Effects 

Due to thermal pressure support, the maximum gravitational mass of a PNS in a CCSN envi- 
roimient can be significantly higher than the cold neutron star (CNS) value [6]. In extreme models 
with soft EOS, or in models with a hyperonic EOS, it is also possible that the baryonic mass of the 
PNS can exceed its CNS value. An interesting case arises when, during a core-collapse event, the 
shock is successfully reenergized but enough material has accreted onto the PNS so that it exceeds 
its maximum CNS baryonic mass. Initially stable, the PNS cools via neutrinos. Since baryonic 
mass is conserved, once sufficient thermal support is lost, no stable configuration exists and the 
PNS will collapse to a BH. We perform simulations with the ul5 model from [12] and we artifi- 
cially adjust the amount of neutrino heating via the scale factor /heat (see [5] for details) to achieve 
explosions. The u75 is an extreme model. In a failing CCSN, it forms a BH roughly ~300 ms after 
bounce with a baryonic (gravitational) mass of ~2.62 (~2.44)M0. For comparison, the LS220 
EOS has a maximum CNS baryonic (gravitational) mass of ~2.41 ('^2.O4)M0. Due to the pres- 
ence of a compositional interface located at a baryonic mass coordinate of ~ l.SMa,, where the 
density drops by ~50%, it is possible to launch explosions very close to the time of PNS collapse 
and achieve the scenario discussed above. In the right panel of Fig. 1, we show the baryonic mass 

^NNBrem is a more dominant contributor to the neutrino luminosity in lower ZAMS mass progenitors, i.e. 
SMq < M < 20 M;,. In these models the temperatures are typically lower and NNBrem, which scales = T^'^, becomes 
comparable to other neutrino pair-production processes, which scale oc T^. 
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of the PNS at BH formation for 7 different values of /heat- The left panel shows the shock position 
until BH formation for these same models. In the /heat = 1.30, 1.32, and 1.325 models, the PNS 
never ceases to accrete material. In the remaining models, as /heat increases, the shock is reener- 
gized at earlier times leaving PNSs with increasingly smaller baryonic masses. As neutrinos extract 
thermal energy from the PNS, the latter becomes unstable and collapses to a BH. The timescale 
for neutrino cooling is typically 10-100 s, but in the case of the extreme ul5 model, BH formation 
happens within ~ 1 s for these values of /heat- For this particular model, we find that NNBrems 
does not have an observable effect due to the extremely high temperatures (50 — 100 MeV) in the 
outer regions of the PNS. As mentioned in §3.1, very high temperatures lead to the dominance of 
other neutrino pair-production processes over NNBrems at early post-explosion times. As an aside, 
hyperonic EOS can have a much lower maximum CNS baryonic mass but still support large PNS 
baryonic masses when the temperature is high [10]. In this case, the above scenario can occur for 
a much larger range of progenitors and extreme models do not need to be invoked. 
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